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Abstract 

By using the recent experimental measurements of B — > tttt and B — > Kir 
branching ratios, we find that factorization is unable to reproduce the observed 
BRs even taking into account the uncertainties of the input parameters. Charm- 
ing and GIM penguins allow to reconcile the theoretical predictions with the data. 
Because of these large effects, we conclude, however, that it is not possible, with 
the present theoretical and experimental accuracy, to determine the CP violation 
angle 7 from these decays. Contrary to factorization, we predict large asymme- 
tries for several of the particle-antiparticle BRs, in particular BR{B + — ► K + tt°), 
BR(B c i — > K + ir~) and BR(Bd — ► ir + ir~). This opens new perspectives for the 
study of CP violation in B systems. 



1 Introduction 



The theoretical understanding of non-leptonic two body B decays is a fundamental 
step for testing flavour physics and CP violation in the Standard Model and for 
detecting signals of new physics Q— 1|. The increasing accuracy of the experi- 
mental measurements at the B factories |], |7| calls for a significant improvement 
of the theoretical predictions. In this respect, important progress has been re- 
cently achieved by systematic studies of factorization made by two independent 
groups H ||]. These studies, while confirming the physical idea []l0| that factor- 
ization holds for hadrons containing heavy quarks, thq ^> Aqcd, give the ex- 
plicit formulae necessary to compute quantitatively the relevant amplitudes at the 
leading order of the AqcdI^q expansion. They also examine some of the con- 
tributions entering at higher order in Aqcd/^q- The question which naturally 
arises is whether in practice the power-suppressed corrections, for which quanti- 
tative estimates are missing to date, may be phenomenologically important for B 
decays. This problem was previously addressed in refs. |Tl], [l^, |u|. In particular, 
the main conclusion of refs. |ll]] was that non-perturbative penguin contractions of 
the leading operators of the effective weak Hamiltonian, Qi and Q2, although for- 
mally of C(Aqc , £)/?tiq), may be important in cases where the factorized amplitudes 
are either colour or Cabibbo suppressed. The most dramatic effect of these non- 
factorizable penguin contractions manifested itself in the very large enhancement 
of the B — > Kir branching ratios, as was also emerging from the first measurements 
by the CLEO Collaboration [14|. In this case, the effect was triggered by Cabibbo- 
enhanced penguin contractions of the operators Q\ and Q2, usually referred to as 
charming penguins. Since the original publications, about three years ago, several 
other decay channels have been measured ]i~5| , 16, 17\ and the precision of the mea- 
surements is constantly improving. With respect to previous analyses, it is now 
possible to attempt a more quantitative study of charming penguin effects and of 
the corrections expected to the factorized predictions. We now present the main 
conclusions of our new analysis. 



Factorization with \V u b\ and 7 from other determinations 

Using the available experimental information on \V u b\ and on the CP angle 7 
provided by the unitarity triangle analysis (UTA) the branching ratios pre- 
dicted with the factorized amplitudes, including the 0(a s ) corrections computed 
according to ref. ||, fail to reproduce the experimental B — > Kir branching 
ratios that are systematically larger than the theoretical predictions. In addi- 
tion, BR(Bd — > 7r + 7r~), which depends on the semileptonic form factor /^(O), 
is about a factor of 2 larger than its experimental value q We note that the 
value of BR{Bd —* ir + ir~) within factorization is essentially fixed by the measured 
BR{B + —* 7r + 7T°) rate. Thus, contrary to the statement of ref. |J, the predicted 
value of BR(Bd — > ir + ir~) is independent of the theoretical assumptions on the 
value of /tt(0). This holds essentially true also for the B — > Kir BRs since the value 

1 Unless explicitly stated the BRs always refer to the average of particles and antiparticles, e.g. 
BR{B d -> K Q ir a ) = \/2{BR{B° d -> K°ir°) + BR{B° d -> K°ir )). 
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of the "semileptonic" form factor at zero momentum transfer fxi®) is correlated 
to /tt(0) by the approximate SU(3) symmetry. 



Factorization fitting 7 

Even if one ignores the value of 7 from UTA, which is only justified if there are 
contributions to AF = 2 mixing due to physics beyond the Standard Model, 
there are serious difficulties in reproducing the experimental results. In particular, 
BR(Bd — ► K °7T ) and BR{B + — ► K°n + ) are much smaller than their experimental 
values. Moreover, the value of 7 extracted from a fit to the data, 7 = (163 db 12)°, 
is in total disagreement with that from the UTA, 7 = (54.8 ± 6.2)° In 
addition, in order to enhance the B — * Kit rates, the preferred values of /jf (0) = 
0.40 ± 0.02 and /^(O) = 0.34 ± 0.01 are incompatible with the latest theoretical 
estimates, f n (0) = 0.26 ± 0.05 ± 0.04, /^(0)// 7r (0) = 1.21 ± 0.09±g;gg @ and 
f„(0) = 0.28±0.05, fK(0)/U(.0) = 1.28lg:io & wh ereas \V ub \ must have a rather 
low value, \ V u b\ = (2.79 ±0.19) x 10~ 3 instead of that extracted from inclusive ^] 
and exclusive ||] semileptonic B decays, \V ub \ = (3.25 ± 0.29 ± 0.55) x 10~ 3 . We 
conclude that, even relaxing the constraint on 7, it is very difficult to reconcile 
the predictions from factorization with the experimental and theoretical findings. 
For this reason any attempt to extract, within factorization, the value of 7 from 
ratios of BRs, for which the discrepancies with the experiments can be accidentally 
hidden, is not very useful. We think that a preliminary step is to understand the 
missing dynamical effects. 



Factorization and charming penguins 

The inclusion of charming penguin effects, which will be explained in detail in 
sec. ||, considerably improves the situation for the B — > Kir channels, with values 
of \V u b\ and 7 well compatible with other determinations. In contrast to the B — ► 
Kn case, charming penguins are not Cabibbo enhanced in B — » irir decays and 
are thus expected a priori to play a minor role. For this reason they should 
be consistently neglected, together with all other Agc^/m,;, corrections. This 
would leave the problem of a too large predicted BR{B^ — > it + it~) unsolved. A 
natural question is then whether the inclusion of ^QCDl^b effects in B& — > ir + ir~ 
can improve the agreement of the predictions with the experimental data. In 
particular, besides the charming penguins, penguin contractions of Qf and Q% 
(GIM penguins in the notation of ref. Ill]]), which are Cabibbo suppressed in 
B — > Kit, might play an important role. We show that, for numerical values 
of the charming and GIM penguin amplitudes of the expected size, h.QCD/ m b ~ 
0.1-0.2, we can easily reproduce the experimental data for both B — * Kir and 
B — > 7T7T decays while respecting the constraints from the UTA. The sensitivity of 
B — ► 7r + 7r~ to A/ nib effects casts serious doubts on the possibility of extracting 
sin 2a from the coefficient of the sinAm^.t term obtained from CP asymmetry 
measurements. On the other hand, we find that the value of the rate asymmetry, 

Am : rr+rr^ - ~ - (U 

^ d ' ~ BR(B° d - 7T+7T-) + BR(B° d - 7T+7T-) ' 1 ' 
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could be unexpectedly large and call our experimental colleagues for separate mea- 
surements of the B and B BRs. In particular, we find \A{B^ — > K^tt )] = 
0.18±0.06, \A(B d -> ^±71^)1 = 0.17±0.06. We also find \A(B d -> vr+7r-)| = 0.30- 
0.50. In the latter discussed in the following, the results are subject to other 

effects on which we do not have control. For this reason we do not quote an error. 
We simply signal that there is room for a large asymmetry also in B d — ► ir + ir~ 
decays. 



2 Results 

In this section we describe and discuss more in detail the different cases which have 
been considered in our analysis. 

The physical amplitudes for B — > Kit and B — > tttt are more conveniently 
written in terms of RG invariant parameters built using the Wick contractions of 



the effective Hamiltonian [23]. In the heavy quark limit, following the approach of 
ref. ||, it is possible to compute these RG invariant parameters using factoriza- 
tion. The formalism has been developed so that it is possible to include also the 
perturbative corrections to order a s , i.e. at the next-to-leading order in pertur- 
bation theory. We present results obtained with this formalism with the addition 
of the non-perturbative \QCD/ m b corrections to factorization described below in 
this section. An alternative framework is provided by the approach of ref. ||. 
This method differs in the treatment of the 0{a s ) terms; unlike the method of 
ref. ||], the calculations are only valid at the leading logarithmic order and it is 
not clear how the independence of the final result from the renormalization scale 
of the operators of the effective Hamiltonian is recovered. Moreover the Sudakov 
suppression of the endpoint region, advocated in ||, is still rather controversial 
from both the theoretical and phenomenological point of view. For these reasons 
we prefer to postpone the analysis with the approach of ref. || until the theoretical 
situation will become clearer. 

In the leading amplitudes, we have taken into account the SU(3) breaking terms 
by using the appropriate decay constants, fx and /„-, and form factors, /if (0) and 
/tt(0). Strictly speaking, the form factors should be evaluated at the invariant mass 
of the emitted meson (/^(m^), fn{rn 2 K ) or / T (m^)). The difference is however of 
higher order in Agc^/m?, and not Cabibbo or colour enhanced and can safely be 
neglected (it is also numerically immaterial) p4}| . As for ^QCD/fnb corrections, 
we have assumed instead SU(3) symmetry and neglected Zweig-suppressed con- 
tributions. In this approximation, by SU(3) symmetry one can show that all the 
Cabibbo-enhanced Agc^/m^ corrections to B — ► Kit decays can be reabsorbed 
in a single parameter P\. Several corrections are contained in P\: this parameter 
includes not only the charming penguin contributions, but also annihilation and 
penguin contractions of penguin operators. It does not include leading emission 
amplitudes of penguin operators {Qz~Q§) which have been explicitly evaluated 
using factorization. Had we included these terms, this contribution would exactly 
correspond to the parameter Pi of ref. [23]. The parameter Pi (Pi) encodes auto- 
matically not only the effect of the annihilation diagrams considered in [|25|| , but 
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all the other contributions of 0(Aqcd / m b) with the same quantum numbers of 
the charming penguins. In this respect it is the most general parameterization of 
all the perturbative and non-perturbative contributions of the operators Q5 and 
Q& (Q3 an d Q4), including the worrying higher-twist infrared divergent contribu- 
tion to annihilation discussed in ref. (2f|. The parameter Pi is of 0{kQCD/fnb) 
and has the same quantum numbers and physical effects as the original charming 



penguins proposed in [11], although it has a more general meaning. In some of the 
previous analyses, see for example |27|, penguin contractions of the operator Qq, 
computed by using perturbation theory and factorization, are enhanced by taking 
a low effective scale for a s . This procedure produces a physical effect similar to 
that coming from the non-perturbative charming penguins that we are using here, 
since they have the same quantum numbers. 

If one also includes B — > tttt decays we have several other parameters, for 



example P^ IM and P3, in the formalism of ref. [23]. A closer look to P3 shows that 
this term is due either to Zweig suppressed annihilation diagrams (called CPA and 
DPA in ref. |ll|]) or to annihilation diagrams which are colour suppressed with 
respect to those entering Pi. For this reason we have put P3 to zero. P^ IM will be 
discussed later on. 

We give now the explicit expression of the B d — > K + ir~ amplitude as an il- 
lustrative example. In terms of the parameters defined in this amplitude 
reads 

A(B d -» K+7T-) = - V us V: b (E 1 (s,u,u;B d ,K+,ir-)-Pr i (s,u;B d ,K + ,iT~)) 

+ V ts V t tP 1 (s,u;B d ,K + ,7r-). (2) 

Using the approach of || , we have 

E 1 (s,u,u;B d ,K + ,TT~) = a?(i^)(Q?) fact + ^(Kvr)(Q 2 ") fact + E x 

6 

P 1 (s,u;B d ,K + ,7r-) = J2 a i( K *)(Qi)f^ +P l 

i=3 

P? 1M (s,u;B d ,K+,ir-) = - aJ , (^))(Q i ) fect + Pf IM , (3) 

i=3 

where (Qi)f ac t denotes the factorized matrix element, and the parameters a, are 
defined in j|]. The tilded parameters represent Agco/mft corrections; in B — » Kir 
channels the only Cabibbo-enhanced correction is given by P\. This term has no 
arguments since we take it in the SU(3) symmetry limit. 

We use input parameters (like p, fj, the form factors) with errors, and extract 
output quantities (like the BRs, the asymmetries, but also 7, or the form factors 
when they are not used as inputs) with their uncertainties. Let us explain how 
we used the input errors and extracted the output uncertainties. We proceed with 
the usual likelihood method, by generating the input quantities weighted by their 
probability density function (p.d.f.). In the case of theoretical quantities this is as- 
sumed to be flat, whereas the experimental quantities are extracted with Gaussian 
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distributions. Probability density functions, averages and standard deviations are 
then obtained by weighting the output quantities by the likelihood factor 



jC = e -iE,(^-^r p ) 2 /-! ) ( 4 ) 

where cij are the standard deviations of the experimental BRs, BR^ xp , given in ta- 
ble 0. In cases where the experimental input has a systematic error dominated by 
theoretical uncertainties, we should extract the latter with a flat distribution [|lS|], 
We have instead combined the errors in quadrature and extracted all the exper- 
imental quantities with gaussian distributions. Within the present accuracy, and 
taking into account the unknown non-perturbative parameters, this procedure is 
fully justified. We have also verified that by extracting the theoretical errors with 
a gaussian distribution, we obtain very similar results. For more details on the 



likelihood procedure, the reader is referred to [18|, where all aspects are discussed 
at length. 



Results with factorization 

We start by considering the case in which we use factorization and take the CKM 
parameters \V u j,\ and 7 from other experimental determinations. We discuss first 
BR(B + —* 7T + 7r°) since in this case, due to isospin symmetry, we do not have the 
complications due to penguin contractions. Thus, at fixed \ V u b\, the prediction for 
BR{B + — ► 7T + 7r°) only depends on /^(O) (trivial dependences as from f w will be 
omitted in this discussion). By using the theoretical estimate and uncertainty of 
/tt(0) from (nj, and taking into account the uncertainties on \V u b\, we predict in 
this case BR(B + — > 7r + 7r°) = (5.0 ± 1.5) x 10 -6 in very good agreement with the 
experimental average given in table |]. A complementary exercise is to use as input 
\V u b\ and the experimental value of BR(B + — ► 7r + 7r°) in order to extract the value 
of /tt(0). In this case we find f n (0) = 0.28 ± 0.06, in very good agreement with 
lattice and QCD sum rules estimates. This exercise shows that we do not need to 
rely on theoretical calculations for the form factors. Indeed also for /if(0) we only 
need /a - (0)7/^(0) which cannot differ too much from one. Moreover it is likely 
that a large part of the uncertainties of the theoretical predictions cancel in this 
ratio. 

Here and in all the other cases where \V u b\ and 7 are taken from other experi- 
mental determinations, we use as equivalent input parameters the values of p and 
f) given in table [l] from the UTA analysis of ref. |l8|]. These values correspond to 

7 = (54.8 ±6.2)°. (5) 

By using f n (0) either from theory or from the fit to BR{B + — * 7r + 7r ) and as- 
suming factorization, we then predict BR(B c i — > ir + ir~) as a function of 7 only. 
Besides, in order to analyze all B — » Kir decays, we only need /a p (0)// 7I -(0) to 
which the previous considerations apply. Alternatively we may take only \V u b\ 
from the experiments and fit the value of 7. In the first case, the results are given 
in table [2| labeled as "7 UTA" and show a generalized disagreement between pre- 
dictions and experimental data. In the second case, the value of 7 is fitted and 
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/*(<>) 


0.27 ±0.08 




1.2 ±0.1 


P 


0.224 ±0.038 


V 


0.317 ±0.040 


BR(B d -> fsr 7r°) 


10.4 ±2.6 


BR(B + -> K+tt ) 


12.1 ± 1.7 


-> if 7T + ) 


17.2 ±2.6 




17.2 ± 1.6 




4.4 ±0.9 


RR(fi+ -> 7T+7T ) 


5.2 ±1.7 



Table 1: Input values used in the numerical analysis. The form factors are taken from 
refs. jl^ |?Z]/ ; i/ie CKM parameters from ref. jJW and the BRs correspond to our average 
of CLEO, BaBar and Belle results fiTh] , \Tb\ \HjJ . All the BRs are given in units of 10~ 6 . 



BR 


7 UTA 


7 free 


BR 


7 UTA 


7 free 


K%° 


5.9 ±0.2 


5.7 ±0.4 




4.8 ±0.2 


9.1 ±0.5 


K°n + 


11.7 ± 0.5 


11.6 ±0.8 


K+TT- 


9.8 ±0.4 


17.7 ± 1.0 


7T + 7T~ 


8.5 ±0.3 


5.1 ±0.7 


7T + 7T 


4.2 ±0.2 


5.4 ±0.6 


7T 7T 


0.19 ±0.01 


0.59 ±0.04 









Table 2: Results for the BRs obtained with factorization without charming or GIM 
penguins. All the BRs are given in units o/10 -6 . 

the results are labeled as "7 free". In this case the disagreement is reduced for 
BR(B + -> K + ir°) and BR(B d -> K + ir~), and also for BR(B d -» vr+Tr"), but it re- 
mains sizable for BR(B d -» i^°7r°) and BR(B + -» J fC°vr+). The pattern BR(B + -» 
K°TT + ):BR{B d -> K + -K-):BR{B d -> K°ir°):BR{B + -> if +7r°)=2:2:l:l ; which is 
suggested by the data, and is well reproduced when the contribution of the charm- 
ing penguins is large, as discussed in the following, is lost in this case. Moreover 
the fitted value of 7 = (163± 12) is in striking disagreement with the results of the 
UTA. Although one may question on the quoted uncertainty of the UTA result, it 
is clearly impossible to reconcile the two numbers. Thus either there is new physics 
or AQCD/mb corrections are important. We now discuss the latter possibility. 

Factorization with Charming and GIM penguins 

We now discuss the effects of charming penguins, parameterized by Pi. P\ is a 
complex amplitude that we fit on the B — > Kir BRs. In order to have a reference 
scale for its size, we introduce a suitable "Bag" parameter, B%, by writing 

Gf 

Pi =77= UU (0)5i B 1 , (6) 

where Gp is the Fermi constant. We use f w (0) for both B — * Kit and B — * tttt 
channels since, as mentioned before, for charming penguins we work in the SU (3) 
limit, gi is a Clebsh-Gordan parameter depending on the final Ktt (tttt) channel. 
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Figure 1: p.d.f. for <fi, in the case where only P± (left) and both Pi and P 1 GIM (right) are 
included. 

In the case where \V u b\ and 7 are taken from the UTA, by fitting the B — > Kir 
channels and B + — ► 7r + 7r° only, we find 

|J3i| = 0.14 ±0.05. (7) 

Note that the size of the charming penguin effects is of the expected magnitude. 
As for the phase eft = Arg(B\), it is very instructive to consider its distribution, 
which is displayed in fig. [I]: the preferred value of <f> has a sign ambiguity since 
we are fitting the average of the B® and B® BRs (or of the B + and B~ BRs). 
The ambiguity can be resolved by measuring separately particle and anti-particle 
BRs. By using the distribution on the left of fig. |lj, we compute the mean value 
of \(p\ with the result \cj>\ = (75 ± 44)° and leave the sign undetermined. This 
is a reasonable procedure, given the approximate symmetry of the distribution 
and the large uncertainty. In view of the discussion of the particle-antiparticle 
asymmetry which we present at the end of this paper, we note here that the value 
of 4> could be rather large. In table || we give the corresponding predicted values 
and uncertainties for the relevant branching ratios (label "Charming"). We observe 
a remarkable improvement for the Kir channels and a large shift in the value of 
BR{Bd — > 7r 7r°) 0, in spite of the fact that in the latter case penguin effects 
are not Cabibbo enhanced (the 7r°7r° amplitude is however colour suppressed). 
The predicted value for BR(Bd — > tt + it~) remains however much larger than the 
experimental one. 

If one fits the B — ► Kir channels, B + — ► 7r + 7r° and B^ — » 7r + 7r~ simultaneously, 
one finds a better agreement for BR(Bd — > tt + it~) but a rather small value for 
BR{B + — > 7T + 7r°) (column "Charming with tt + tt~" of table |3|). This happens at 
the price of reducing the fitted value of the form factor, /^(O) ~ 0.22, which is 
pushed down by BR(Bd — * 7r + 7r~). In fact the latter has an experimental error 
much smaller than BR(B + — > ir + ir°), and therefore governs the fit. However, we 
do not think that this is the correct procedure: theoretically, BR(B + — ► 7r + 7r°) 
is on much more solid grounds than BR(B c i — > tt + tt~), since it is not affected by 
penguins or annihilations, and thus is much more suitable to constrain /^(O). 

2 This effect was already noticed in jflj . 
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Figure 2: p.d.f. for the largest CP asymmetries, in the case where only Pi (left) and 
both Pi and P^ IM (right) are included. From top to bottom, we give A(B + — > K + tt°), 
A(B d -> K+tt-) and A(B d -> tt+tt-). 



BR 


Charming 


Charming 


Charming 


BR 


Charming 


Charming 


Charming 






With 7T + 7T~ 


+ GIM 






With 7T + 7T~ 


+ GIM 


K°n° 


9.2 ± 1.1 


8.7 ±0.9 


8.8 ± 1.0 


K + 7l° 


9.2 ±0.8 


9.3 ±0.7 


9.3 ±0.7 


K°7l + 


18.3 ±2.1 


17.4 ± 1.8 


17.6 ± 1.8 


K+iT- 


18.2 ± 1.4 


18.6 ±1.4 


18.4 ± 1.3 


7T + 7T _ 


9.1 ±2.5 


5.1 ± 1.8 


4.7 ±0.8 


71 + 71° 


4.8 ± 1.4 


2.7 ±0.5 


3.5 ±0.9 


7T 7T 


0.37 ±0.05 


0.36 ±0.05 


0.69 ±0.30 











Table 3: BRs with charming or charming and GIM penguins. All the BRs are given in 
units of 10~ 6 . 
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In order to reduce the predicted BR(Bd — ► vr + vr~) without affecting BR(B + — ► 
7r + 7r°), one may include other effects of the same order of the charming penguins, 
as for example the GIM penguins introduced in ref. [11]. In this case we fit all the 
BRs given in table |l]. With GIM and charming penguins included, we find 

|2?i| = 0.16 ±0.03, |0| = (56 =b 32)° , 
\Bf m \ = 0.23 ± 0.11 , |0 GIM | = (135 ± 37)° , (8) 

where the notation is self-explaining. We have given the absolute value of since, 
as in the previous case, the sign ambiguity persists when we include GIM pen- 
guins. The distribution is also shown in fig. [|. The results for the BRs can be 
found in table || with the label "Charming+GIM" . They show that the extra GIM 
parameter improves the agreement for the measured B — * tttt BRs. We do not 
claim, however, to be able to predict BR{Bd — > vr + vr~): our results instead show 
that accurate predictions for B& — > vrvr decays can only be obtained by controlling 
quantitatively the 0(AQCD/ r nb) corrections, which is presently beyond the theo- 
retical reach. Estimates for charming penguin effects can also be obtained by using 
some phenomenological model, as for example done in ref. fllfl . We observe that 
the sensitivity of the BRs to the value of 7 is lost, with the present experimental 
accuracy, once penguin effects are introduced. Indeed when one tries to fit B\ 
(B± IM ) and 7 simultaneously, one finds that the value of 7 is essentially undeter- 
mined. From the above discussion it clearly emerges that one of the important 
step for the improvement of this kind of analyses is a more precise measurement 
of BR{B+ -> vr+vr°). 



Particle— Antiparticle asymmetries for the Branching Ratios 

The large absolute values of <p, and the sizable effects that penguins have on 
the BRs, stimulated us to consider whether we could find observable particle- 
antiparticle asymmetries as the one defined in eq. (|l|). We find large effects in 
BR(B+ -> K+7T ), BR{B d -» K + ir~) and BR(B d -> vr+vr-), as shown in fig. g 
As discussed before, for BR(Bd — ► vr + vr~) our predictions suffer from very large 
uncertainties due to contributions which cannot be fixed theoretically. For this 
reason, the values of the asymmetry reported in table ||] are only an indication that 
a large asymmetry could be observed also in this channel. The sign ambiguity of 
<j> is reflected in the asymmetry A ~ sin 7 sin 0. This ambiguity can be solved only 
by an experimental measurement or, but this is extremely remote, by a theoreti- 
cal calculation of the relevant amplitudes. For each channel, we give the absolute 
value of the asymmetry in table |I[ Note that within factorization all asymmetries 
would be unobservably small, since the strong phase is a perturbative effect of 
0(a s ) M. The possibility of observing large asymmetries in these decays opens 
new perspectives. These points will be the subject of a future study. 
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\A\ 


Charming 


Charming ± GIM 


\A\ 


Charming 


Charming ± GIM 


K°7l° 


0.02 ±0.01 


0.05 ±0.03 


K + 7l° 


0.23 ±0.10 


0.18 ±0.06 


K°n+ 


0.00 ±0.00 


0.03 ±0.03 


K+ir- 


0.21 ±0.10 


0.17 ±0.06 


7T + 7T~ 


0.36 ±0.16 


0.52 ±0.18 




0.40 ±0.19 


0.58 ±0.29 



Table 4: Absolute values of the rate CP asymmetries for B — > Kit and B — > ixix decays. 
The columns labeled by "Charming" and "Charming + GIM" correspond respectively to 
the cases in which only Pi and both Pi and are introduced. The asymmetry in 

B — > 7r + 7T° vanishes exactly. 

Conclusion 

We have analyzed the predictions of factorization for B — > rnr and B — > Kir 
decays. We note that the normalization of all the other BRs is essentially fixed 
by the value of BR(B + — > 7r + 7r°) and SU(3) symmetry. Even taking into account 
the uncertainties of the input parameters, we find that factorization is unable to 
reproduce the observed BRs. The introduction of charming and GIM penguins []ll]] 
allows to reconcile the theoretical predictions with the data. It also shows however 
that it is not possible, with the present theoretical and experimental accuracy, to 
determine the CP violation angle 7. Contrary to factorization, we predict large 
asymmetries for several of the particle-antiparticle BRs, in particular BR(B + — > 
K + 7t°), BR(Bd — > K + 7r~) and BR(Bd — > ir + ir~). This opens new perspectives 
for the study of CP violation in B systems. 
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